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ABSTRACT 

Aims. This study of SDSS J080434.20+5 10349.2 is primarily concerned with the double-hump shape in the light curve and its 
connection with the accretion disk in this bounce-back system. 

Methods. Time-resolved photometric and spectroscopic observations were obtained to analyze the behavior of the system between 
superoutbursts. A qualitative geometric model of a binary system containing a disk with two outer annuli spiral density waves was 
applied to explain the light curve and the Doppler tomography. 

Results. Observations were carried out during 2008-2009, after the object's magnitude decreased to V ~ 17.7 ±0.1 from the March 
2006 eruption. The light curve clearly shows a sinusoid-like variability with a 0.07 mag amplitude and a 42.48 min periodicity, which 
is half of the orbital period of the system. In September 2010, the system underwent yet another superoutburst and returned to its 
quiescent level by the beginning of 2012. This light curve once again showed a double-hump-shape, but with a significantly smaller 
(~ O.Olmag) amplitude. Other types of variability like a "mini-outburst" or SDSS1238-like features were not detected. Doppler 
tomograms, obtained from spectroscopic data during the same period of time, show a large accretion disk with uneven brightness, 
implying the presence of spiral waves. 

Conclusions. We constructed a geometric model of a bounce-back system containing two spiral density waves in the outer annuli of 
the disk to reproduce the observed light curves. The Doppler tomograms and the double-hump-shape light curves in quiescence can 
be explained by a model system containing a massive >0.7Mq white dwarf with a surface temperature of ~12000K, a late-type brown 
dwarf, and an accretion disk with two outer annuli spirals. According to this model, the accretion disk should be large, extending to 
the 2:1 resonance radius, and cool (~2500K). The inner parts of the disk should be optically thin in the continuum or totally void. 

Key words, stars: - cataclysmic variables - dwarf nova, individual: - stars: SDSS J080434.20+5 10349.2 - accretion, accretion disks 
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1. Introduction 

A widely accepted evol utionary theory of catac lysmic variables 
(CVs), as presented in iKolb & BaraffS (fT999l and references 
therein), predicts a significa nt accumulation o f systems around 
the orbital period minimum ('Paczynski"l98lV It also envisions 
that ~70% of the cuiTent CV population has evolved past the 
orbital period minimum and formed so-called bounce -back sys- 
tems that contain a brown dwarf as a secondary (iPatterson I 
1201 lb . According to the calculations, the orbital period of a 
bounce-back system should be slightly longer than before reach- 
ing the period minimum. Bounce-back systems are expected 
to float within the 80 - 90 min orbital period range. However, 
neither a peak in the period distribution nor any past period- 
minimum systems were observed prior to the Sloan Digital Sky 
Serve (SDSS). The SDSS has identified more than a hundred 
new faint CV systems w ith accurately measured orbital peri- 
ods dGansicke et al.ll2009l) . A significant portion of the new CVs 
have periods clustered at the lower end of the orbital period 
distribution and some of them are good bounce-back system 
candidates. SDSS J103533.03-H0 55158.4 (hereafter SDSS 1035) 
was proposed by iLittlefair et al.l ( 12008.) as a solid example of a 
bounce-back system based upon the system parameters found 
using aj)arametrized model of an eclipse observed in its light 
curve.lAviles et al. (2010) included SDSS J123813.73-033933.0 
(hereafter SDSS 1238) in the list of bounce-back CVs using a di- 



rect IR broad band photometry of the secondary and modeling 
the spectral energy distribution of the system. 

The subject of this paper is SDSS J080434.20+5 10349.2 
hereafter SD SS0804), which has a very peculiar light curve 
ISzkodv et al.l i2006). A similar light curve has only been ob- 
served in one oflier system of SDSS 1238. Hence, SDSS0804 is 
considered a twin to SDSS1238. Since its discovery, SDSS0804 
underwent tw o superoutbursts , confi rming its membership as a 
WZ Sge-type. IZharikov etd] (l2008h included it in the list of 
bounce-back candidates based upon its orbital period and mass 
ratio estimates. SDSS0804 exhibits some phenomena ("bright- 
enings" and "mini-outbursts") that are not common in other 
CVs ( Zharikov et al. 2008). Its superoutburst in September 2010 
(iPavlenko et al.l201 ll) . which was only about four years after the 
previous superoutburst, was unexpected. A general description 
of the characteristi cs and observation al features of SD SS0804 
can be found in Szkodv et al. (2006), Pavlenko et al 



IZharikov et al. 
(1201 Ih . 



(200 8), K ato et al. 



(200' 

(2009), and Paylenko et^ 
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In this paper we report the results of monitoring SDSS0804 
in 2008-2009 and 2012, during periods when the object almost 
retuned to its pre-superoutburst, quiescent level. The main aim 
of this study is to investigate the system's behavior to determine 
the structure of the accretion disk during the quiescence. 

In Sect. |2] we describe our observations and the data reduc- 
tion. The light curve analysis and results of infrared and spectral 
observations are presented in Sects. [3]and|4l A discussion of the 
accretion disk structure of this bounce-back system is given in 



1 



Zharikov, S., et al.: The accretion disk in the post period-minimum CV SDSS0804 



T 1 1 r 



"1 1 1 r 



T 1 1 r 



"1 1 1 r 



"1 1 1 r 



12.6min variabiUty 
I > 



I 



SDSS1238-type WZ-Sge -type "mini-outbursts" 
variability superoutburst double-humped 

March 4, 2006 light curve 



J I L 



J I L 



J I L 



WZ-Sge -type 
superoutburst 
Sept. 22, 2010 

J I I I I I I I L 



53500 



54000 



54500 55000 
HJD +2400000 



55500 



56000 



Fig. 1. Light curves of SDSS0804 are shown from 2005 t o 2012 in the V band. The plot is based on a combination of data obtained 
by us, including those published in lZharikov et al.l ( l2008l) . and data accumulated by the AAVSO. 



Table 1. Log of the time-resolved photometric observations of 
SDSSJ080434.20-H5 10349.2. 



Table 2. Log of the time-resolved spectroscopic observations of 
SDSS J080434.20H-5 10349.2. 



Date 


HJD start -l- 


Exp. time Duration 




2454000 


num. of integrations 



Photometry 








2008-02-06 


502.634 


30sx821 


6.8h 


2008-02-07 


503.626 


30sx297 


2.5h 


2008-02-08 


504.608 


30sx416 


3.5h 


2008-03-08 


534.656 


60sx300 


5. Oh 


2008-03-09 


535.643 


60sx225 


3.8h 


2008-03-10 


536.632 


60sx250 


4.2h 


2008-03-11 


537.620 


60sx290 


4.8h 


2008-11-22 


593.511 


120sx90 


3. Oh 


2008-11-23 


594.469 


120sx41 


1.4h 


2009-01-23 


855.724 


90sxl05 


2.6h 


2009-01-24 


856.690 


90sxl78 


4.5h 


2009-01-25 


587.688 


90sxl52 


3.8h 


2009-02-27 


890.705 


60sxll0 


1.8h 


2012-01-13 


1940.786 


60sx 107 


5.6h 


2012-01-17 


1944.884 


60sx 87 


3.1h 


2012-01-19 


1946.881 


60sx 104 


3.6h 


2012-01-21 


1948.931 


60sx 75 


2.4h 


2012-01-22 


1949.781 


60sx 149 


5.5h 


2012-01-25 


1952.780 


60sx 123 


5.5h 


2012-01-26 


1953.711 


60sx 179 


6.7h 


2012-01-27 


1954.906 


60sx 113 


2.6h 


2012-01-28 


1955.836 


60sx 122 


4.3h 


2012-01-29 


1956.835 


60sx 123 


4.3h 


2012-01-30 


1957.855 


60sx 110 


3.8h 



Sects. |6] and |2l The light curve simulation of bounce-back sys- 
tems is described in Sect. [8] and our conclusions are given in 
Sect.|9] 

2. Observations and data reduction 

Time-resolved photometry of SDSS 0804 was obtained using the 
direct CCD image mode at the 1 .5m and 0.84m telescopes of the 
Observatorio Astronomico Nacional at San Pedro Martir (OAN 
SPlvQ) in Mexico. The log of photometric observations is given 



Date 


HJD start + 


Exp. time 


Duration 




2454000 


num. of integrations 




Spectroscopy 








2008-02-06 


502.634 


490sx41 


8.3h 


2008-02-07 


503.626 


490sx53 


8.0h 


2008-03-10 


535.658 


490sx40 


5.8h 


2008-03-11 


536.648 


490sx39 


5.9h 


2008-03-12 


537.641 


490sx40 


6.1h 


2009-01-24 


855.738 


483sx32 


5.3h 


2009-01-25 


856.674 


483sx40 


6.3h 


2009-01-26 


857.630 


483sx40 


6.1h 


2009-01-27 


858.646 


483sx46 


6.9h 


2012-01-25 


1952.670 


500sx52 


7.9h 


2012-01-26 


1953.641 


500sx57 


8.0h 



http://www.astrossp.unam.mx 



in Table [T] Photometric data were calibrated using Landolt's 
standard stars. The errors of photometry were calculated from 
the dispersion of the magnitude of the comparison stars in the 
object field. The errors ranged from 0.01 to 0.05 mag. 

Long-slit observations were obtained with the Boiler & 
Chivens spectrograph and the 2.1-m telescope of OAN SPM. 
The spectra span from 4000 to 7100 A with a resolution of 
3.03 A pixel"' . To improve the signal-to-noise ratio, we obtained 
a series of phase-locked spectra: ten spectra were taken at equal 
phase intervals over a single orbital period (Porb = 80.5 min). 
This sequence of spectra was repeated at exactly the same phase 
intervals for subsequent periods and on subsequent nights. This 
allowed us to calculate the phase-averaged spectra by summing 
up spectra at the same orbital phase obtained during several cy- 
cles without decreasing the time resolution. The log of the spec- 
tral observations is given in Table |2] There is a small difference 
in the exposure times at different epochs because different CCDs 
were used, and each CCD has its own readout time. 

The J band photometry of SDSS0804 was obtained in 
October 2009 with the "CAMILA" IR CCD camera attached to 
the 0.84m telescope of OAN SPM. The data were reduced by the 
standard way and were calibrated using IR standards observed 
before and after the object. 
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Fig. 2. Power spectra (top panel) of SDSS0804 is shown for different observational runs; the corresponding light curves are folded 
with the orbital period (upper axis of the middle panels). The half orbital period (bottom axis of the middle panels) and the 12.6 min 
period (lower panels) are also shown. 
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Fig. 3. Power spectra (top panel) of SDSS0804 using the 2012 
data and the corresponding light curves folded with the orbital 
period (upper axis of the lower panel) and the half orbital period 
(bottom axis of the lower panel) are shown. 



3. Light curve variability 

The light curve of SDSS0804 during 2005-2012 is shown in 
Figure[T] Before the March 2006 superoutburst the brightness of 
the system was V ~ 17.9 and the objec t displa y ed SPSS 1238- 
like long-term variability (S zkodv et al.l (l2006l) : IZharikov et al.l 
(2006, 2008)). On March 4, 2006, a superoutburst commenced, 
exhibiting all attributes of WZ Sge-typ e: a large arnphtud e 
> 6 mag, superhumps, and 11 echoes tPavlenko et al.l 120061) . 
The maximum brightness during the superoutburst reached 
abo ut V ~ 12 m agnitude. Based upon the superhump period 
( Ka to et al.ll2009h. the mass ratio o f the system was estimated to 
be ^ * 0.05 (IZharikov et al.ll2008l) . which is a typical value for 
bounce-back systems. 

About eight months after the first superoutburst, the bright- 
ness decreased to 17.1 magnitude. Variability at different time- 
scales was detected: a rapid modulation with a period of Pi = 
12.6 min and amplitude of ~ 0.05 mag was interpreted as a 
strong nonradial pulsation of the cooling (i.e., after the outburst) 
white dwarf, that entered the instability strip (Pavlenko et al] 
I2OO6); significant short-term variations with periods P2 — 21.7 
min, P3 - 14.1 min, and P4 - 4.28 min that may be ad- 
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Fig. 5. Flux-calibrated time-averaged spectra of SDSS0804 obtained in different epochs. 

The photometric data were analyzed for periodicities us - 
ing a discrete Fourier transformation code dDeemind Il975h . 
The power spectra for different observational runs are given in 
Figures |2] and |3] The power spectra clearly show double or- 
bital frequency in all observational runs. The amplitude of the 
double-hump-shape light curve decreases with brightness (from 
~ 0.07 mag at HJD 24541 10, ~ 0.02 mag at HJD 245456 to only 
~ 0.01 at HJD 2455953). The 12.6 min period, attributed to WD 
pulsations, was only detected in the February - March 2008 and 
January 2009 runs. This period disappeared afterwards (Figs|2] 
andO. The light curves folded by the 42.5 min and 12.6 min pe- 
riods are shown along with their corresponding power spectra in 
Figures|2] and |3] . The long-term SDSS1238-like variability and 
mini-outburst phenomena were not detected during our moni- 
toring of the system between superoutbursts and after the last 
superoutburst (HJD 2454500 — HJD 2456000). 




Fig. 4. J band image obtained in October 2009 using Camila at 
the 0.84m telescope of the OAN SPM. 



ditional pulsation modes (iPavlenko et al.ll20Tll) : and the super- 
humps were replaced with orbital double-hump-sha pe light vari- 
ations dZharikov et al.l 120081: [Pavlenko et al.ll2007h . In addition 
to these, the object also showed a new phenomenon - recurring 
(~ 32 day) "mini-outbursts" that lasted approxi mately four days 
and in creased the brightness to up to ~ 0.6 mag (IZharikov et al.l 
I2008h . 

Given the diversity of the variability displayed by 
SDSS0804, we continued monitoring the system with an aim 
to find any systematics in its behavior In the three years af- 
ter the 2006 superoutburst, the V-band magnitude of the ob- 
ject decreased to about 17.7 mag. However, the object never 
reached the pre-superoutburst level of March 2006 because in 
Septemb er 2010 S DSS0804 suddenly underwent another super- 
outburst (iPavlenko et al.. .20 1 1 ). only ~ 4.5 years after the pre- 
vious one. This short recess between superoutbursts is atypical 
for WZ Sge-type systems. The amplitude of this superoutburst 
was similar to the previous, but the shape of the light curve 
was different: the brightness decrea sed faster and fewer echoes 
were detected dPavlenko et al.ll20Tll) . By December 201 1, more 
than a year after this superoutburst, the brightness of the system 
fell practically to the quiescence level. From December 201 1 to 
January 2012 we obtained a new set of photometric and spectro- 
scopic data at quiescence to compare with our previous observa- 
tions. 



4. Infrared data 

i Kato et al.1 d2009h performed infrared photometry of SDSS0804 
using OAO/ISLE in March 2007, about one year after the super- 
outburst. The IR magnitudes were estimated to be 7 = 17.29(5), 
H - 16.97(5), and K - 16.41(6). During these observations 
the object was ~ 17.1 in the visual band, about 0.7-0.8 mag 
above the pre-outburst quiescence level. Our IR observations 
were made in the J band in October, 2009. The object magni- 
tudes was J - 17.3(1) (see Fig.lUl. Hence, the object did not 
change significantly in the J band although it was still cool- 
ing and becoming fainter in the optical. We assume that the 
bright ness of the object in the IR, as measured by Kato et al] 
d2009h . corresponds to the quiescence state. On timescales com- 
parable to the orbital period, the system is expected to show 
some variability in the near-IR due to the elliptical shape of the 
Roche-lobe-filling secondary. The ellipsoidal variability of the 
secondary can be calculated and it has been taken into account 
as an additional error of about 0. 15 mag in the following consid- 
erations of spectral energy distribution. 

5. Spectral energy distribution 

The optical spectrum of SDSS0804 shows a blue continuum 
with broad absorption lines (caused by the white dwarf) that 
surrounds the double-peaked Balmer emission lines (caused 



4 



Zharikov, S., et al.: The accretion disk in the post period-minimum CV SDSS0804 



-15 



— I 1 1 1 1 1 r- 



I ' ' ' ' I ' ' ' ' I 

r=-i.93 



-1 — I — I — I — I — I — I — I — I — 1— 

Model 

L5 

Accretion disk 

WD I2000K M,„„=0.7 

• JHK 

SDSS 0804 



-15.5 
-16 
-16.5 



^ -17 

60 
O 



-17.5 



-18 



_l I I I I I L_ 



_l I I I I L_ 



3.6 



3.7 



1< 

'e 

o 

M 
u 
u 

Ph 
60 
O 



-15 
-15.5 

-16 - 
-16.5 

-17 
-17.5 



— I 1 1 1 1 1 r- 



3.8 



3.9 



4.1 



4.2 



4.3 



4.4 



log X (A) 



I ' ' ' ' I ' ' ' ' I 
rQ=-2.33 



-1 — I — I — I — I — I — I — I — I — 1— 

Model 

L4 

Accretion disk 

WD 1 2000K M,„„= 1 .0 

• JHK 

SDSS 0804 




.jg I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

3.6 3.7 3.8 3.9 4 4.1 4.2 4.3 4.4 

log X (A) 

Fig. 7. Spectral energy distribution of SDSS0804 is plotted with the components of the best-fit models. The tiny black line is the 
observed quiescent optical spectrum, red dots are JHK IR measurements. Overplotted with a thick black line is the sum of all 
components that comprise the close binary. The green line is the WD spectral model, the blue line corresponds to the accretion disk 
contribution, and the maroon line represents the brown dwarf spectrum. The top figure shows the solution coiTesponding the 
minimum in which the accretion disk dominates in the IR. The bottom figure presents the solution with the standard spectral index 
of r = Fo = -7/3 = -2.33 and the secondary dominating in the IR. 



by the accretion disk). The spectrum undergoes some changes 
throughout the outburst cycle. The flux-calibrated spectra of 
SDSS0804, which cover a period of more than two years prior 
to the 2006 superoutburst until one year after the 2010 super- 
outburst, are shown in the upper panel of Figure|5] As already 
mentioned above, during the 2008-2009 the object was about 
0.3 mag brighter than before the superoutburst. But in January 
2012, the object brightness returned to the pre-outburst level 
of March 2006. The continuum of the post-outburst spectra is 
bluer, Balmer absorption lines are deeper, and emission lines are 
weaker than in the pre-superoutburst spectrum. The equivalent 
width of Hq,, for example, changed from w -90 to « -15 in 2004 
and 2008. In January 2012, the equivalent width of the H„ line 
increased to up to « -70, about four times larger than values 
obtained between superoutbursts. However, the value does not 
reach the level seen in 2004 (see Figure|6]l because the object 
never reached the same quiescent state prior to the first super- 
outburst. 



We fit the spectral energy distribution (SED) o f SDSS0804 
in quie scence with the model described in detail bv lAviles et al.l 
(120101) using a combination of the SDSS spectrum obtained 
in 2004 wi fl i the IR data of Kato et al. (2009). Following 
lAviles et afl (l2010i) . we assumed that the total flux, F*(/l), is 
the sum of the radiation from a hydrogen WD (DA type) 
FwoiTes, A), an accretion disk with Fad ~ and a red/brown 
dwarf with Fbd('I): 

F\A) = FwD(reff, A) + FadW + F^^^iA). (1) 

Brown dwarf fluxes were taken from the literature 
(iMcLean et all I2003L l2007h and online sourcefl The white 
dwarf radii were calculated using the mass-radius relation 

^ see |http://web.mit.edu/ajb/www/browndwarfs/| 
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Fig, 6. Hq, line profile of the object at different epochs is shown. 
The continuum is subtracted. 
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Fig. 8. Relations of distance, white dwarf temperature, spectral 
type of the secondary, mass of the primary, contribution of the 
accretion disk, and;^'^ (from the top to the bottom) vs. the slope 
of the power law spectrum of the accretion disk. The two models 
presented in Figure |7] are marked by vertical lines. 



from iNauenberj (Il972h for an Mwd = 0.6 - 1 . 1 Mq range, 
with steps of a 0. 1 Mq. Spectra of WDs with pure hydrogen at- 
mosphere in the 4000-25000A spectral interval w ere obtained 
using ATLAS9 dKuruczl [1991 and SYNTH (Pis kunovlll992h 
codes for the appropriate range of temperatures. The calculations 
were performed with a 1 000 K step in a temperature range from 
Teff =10000 to 18 000K and with a surface gravity g = y^. 

The spectra were normalized to Aq = 5500A and the contribution 
of the WD is 

FwD(T,ff, A) = CdS) * F^'^(T,n; A), 

where Ci(6) = io-"-4*(v+'5+<)^ V = 17.9 is the object's bright- 
ness in quiescence, and 5 is a parameter that determines the con- 
tribution of the flux from the WD to the total flux in the V band. 
Finally, My - 21.109 is a constant to convert magnitude into 



flux (in ergs/cm^ /s/ A) in the V band. The spectra of the accre- 
tion disk was assumed to be a simple power law 



Fad(/1) = (Ci(0)-Ci(5))x 



^0 



where (Ci(0) - Ci(6)) determines the contribution from the ac- 
cretion disk in the V band, assuming that the WD and the accre- 
tion disk are the only sources in that wavelength because the only 
other contributor is the brown dwarf, which has an only negli- 
gible flux in the V band. F is the slope of the accretion disk's 
power-law spectrum. 

The distance to the object is estimated to be 



d — Ry 



Fbb(reff,5500A) 

' Fwd(5500A) 



where F^^iT^ff, 5500A) is the black body flux at ^ = 5500A 
with an effective temperature Tgff. Observed SEDs of red/brown 
dwarfs with spectral types between M6 to L5, normalized to fit 
the observed flux in the J band, were used. The bolometric cor- 
rection to the J mag nitude for each spectral type was taken from 
iTinnev et alJ (l2003h . 

The free parameters of this three-component model are the 
white dwarf effective temperature Teff , the mass of the white 
dwarf Mwd, the spectral type of the secondary star SpT, the ratio 
6 of WD-to-accretion disk contribution in the V band, and the 
slope of the accretion disk spectrum F. 

Results of the fitting are presented in Figures|2] and [8] A 
marked minimum in occurs at F = -1.93 (see the vertical 
dash-dotted line in Figure[8]l. This solution implies an Mwd = 
0.7 + 0.1 Mo white dwarf with Two = 12000 + lOOOK tempera- 
ture, a ~ L5+2 brown dwarf, and an accretion disk with a 50% 
contribution in the optical spectrum that dominates in the in- 
frared. The system will be located at a distance of about ~ 200pc. 
Indeed, this solution works for a range of -2.2 < F < -1.8 
where the variation is weak . 

If constraints on are relaxed even more, another solution 
becomes viable for the spectral slope of -2.75 < F < -2.2 
and -1.8 < F < -1.3, which inclu des the standard sp ectral 
index of F = Fq = -7/3 = -2.33 (lLvnden-Bel]|ll9690 (see 
the dashed line in Figure[8]l. The fit parameters for solutions in 
these F ranges correspond to a system formed by a more massive 
MwD = 1.0 + 0.1 Mq white dwarf with the same temperature, a 
brown dwarf of slightly earlier spectral type ~ L4^^, and an ac- 
cretion disk whose contribution is only ~ 30+ 10% in the optical 
spectrum dominated by the white dwarf. In the infrared the con- 
tribution of the accretion disk is comparable to the secondary 
star. Accordingly, the system would be located at a distance of 
about 120- 130 pc. 

The difference between these two solutions arises from the 
different contribution of the accretion disk in the infrared. In the 
model with the lowest x^ the spectrum of the accretion disk is 
flatter compared to the standard model, as expect ed for an opti- 
cally thin disk in the continuum (lldan et al ]l2008h . It also means 
more contribution from the cooler outskirts of the disk in the in- 
frared. Therefore, if the disk dominates in the IR, the distance 
estimate to the system increases, which leads to a less massive 
white dwarf with a larger radius and higher luminosity. However, 
there is evidence that sh ort period CVs tend to contain massive 
(~O.9M0) white dwarfs (Littlef air et ani2008h . It is difficult to 
choose the right solution based on quality of the available data 
and the uncertainties in applied models, but it is beyond doubt 
that the contribution of the accretion disk in the continuum of 
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Fig. 9. Trailed spectra around the Hq, line folded with the orbital period of the system and the corresponding Doppler maps. The 
center of mass, both stellar components, the stream of mass transfer, and the circle corresponding to the 2: 1 resonance radius of the 
accretion disk are marked on the tomograms. 



bounce-back systems is far less than in ordinary CVs. Additional 
ultraviolet and infrared spectroscopy would be necessary to de- 
termine fractions of fluxes from the stellar components and sep- 
arate the spectrum of the accretion disk. 

Additionally, we calculated models for spectra obtained im- 
mediately after the superoutburst. The infrared values of the sys- 
tem were kept the same (within the error range) as in quies- 
cence. The temperature inferred for the best-fit model is Twd - 
15000 + lOOOK. This value is not surprising because an increase 
of white dwarf surface temperature during supe routburst has 
been observed recently in a number of other CVs dGodon et al.l 
l2004tlSion et"an[T99Q) . In the case of SDSS0804 the appear- 
ance of WD pulsations after the superoutburst and their sub- 
sequent disappearance is also good evidence of temperature 
change (.Pavlenkoil2009h . 

6. Doppler tomography 



, Jan. 2012 



We used Doppler tomography (iMarsh & Horne' '1988') to study 
the structure of the accretion disk of SDSS0804. We generated 
Doppler maps from our time-resolved spectra of the obj ect us- 
ing the maximum entropy method as implemented by Spruit! 
( Il99 8fl Trailed spectra around the H„ line and their correspond- 
ing Doppler maps are shown in Figure|9] The orbital period of 
Porb - 0.059 days, a white dwarf mass of Mwd = 1 .0 M©, and 
mass ratio of ^ = 0.05 are used to overlay positions of the stel- 
lar components and the stream on the Doppler mapfl The in- 
clination angle / = 70° is arbitrarily chosen as a compromise 

^ http://www.mpa-garching.mpg.de/~henk/pub/dopmap/ 



Using Mw 



0.7 Mq only influences overlaid graphics, not the 
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Fig. 10. Same as shown in Figure 7 but for the Hp line and only 
the 2012 data. 



between strong double-peaked emission lines in the spectra and 
absence of obvious eclipses in the light curve. A possible pres- 
ence of weak ecl ipses during the superoutburst, as reported by 
iKato et al.l (|2009|) . supports the idea of a high-inclination sys- 
tem. Trailed spectra folded with the orbital period of the system 
clearly show an S-wave emanating from the bright spot created 
by the impact of the stream with the accretion disk. The Doppler 
maps show not only the bright spot, but also two distinctive 
bright regions on the ring, indicative of an accretion disk. The 
bright regions repeat from epoch to epoch. The brighter one is 
located, as expected, where the mass transfer stream hits the ac- 



maps; it slightly shifts the position of the Roche lobe and the accretion 



stream. The decrease of the circle radius denoting the 2: 1 resonance is 
insignificant (from ~680 to ~610km/s). 
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Fig. 11. Plots of the mass ratio and the mass transfer rate vs the orbital period of CVs. In the schema (left panel), locations of 
different CV types are marked along with the 3:1 and 2:1 resonance radii. The lower box in that panel corresponds to the location 
of bounce-back systems, as labeled. In the right panel observed systems with known parameters are plotted with the mass ratio axes 
marked on the left and estimated mass transfer rates on the right axes. The bounce-back candidates are enclosed in a dash-dotted 
box. 



Table 3. Parameters of bounce-back candidates 



NN/Object 


Pot* 


V 


q 


M, 


M2 


eff 


i 


LC' 


References 




(days) 


(mag) 




(Mo) 


(Mo) 


(K) 


n 






1. GWLib* 


0.0533 


19.1 


0.060 


0.84 


0.05 




11 




van Soaandonk et al. (2010) 


2. SDSS1238 


0.056 


17.8 




1.0 


«0.09 


12000 


>70 


+q 


Zharikov et al. 2006 


3. V455 And* 


0.0563 


16.5 


0.060 




>M9 


11500 


83 


+q 


Araujo-Betancor et al. 2005 


4. AL Com* 


0.0567 


20.0 


0.060 






16300 




+q 


Patterson et al. fl996) 


5. SDSS1035 


0.057 


18.7 


0.055 


0.94 


0.05 


10100 


83 




Littlefair et al. (2008) 


6. SDSS0804* 


0.059 


17.9 


0.05 


1.0(1) 


0.05 


12000 


>70 


+q 


Zharikov et al. 2008 


7. EG Cnc* 


0.060 


18.8 


0.035 






12300 




+s 


Patterson et al. (1998) 


8. RX1050-14 


0.062 


17.6 


<0.055 






13000 


<65 




Mennickent et al. (2001) 


9. GD552 


0.0713 


16.6 


<0.052 




<0.08 


10900 


<60 




Unda-Sanzana et al. (2008) 


10. RE1255 


0.083 


19.0 


<0.06 


>0.9 


<0.08 


12000 


< 5 




Patterson et al. (2005J 


11. SDSS1610** 


0.0582 


19.0 












+q 


Mukadametal. (2010) 



' light curve (LC) features: "+" LC shows a double-hump during the orbital period; 
"s" - during superoutburst; "q" - during quiescence; "-" absence of double-humps in LC. 
* - objects that demonstrate WZ Sge-type superoutburst. 

** The mass ratio of SDSS1610 is not known, but the object shows similar observational characteristics to selected candidates and the double- 
humped light curve. 



cretion disk. This bright region is commonly seen in most CVs, 
but here it looks extended along the mass transfer stream line. 
The other bright region is extended and is located completely on 
the opposite side of the disk atV ~ 750km/s. This bright region 
is not common for accretion disk systems bu t is a feature shared 
by the bounce-back candidates SDSS1238 (lAviles et al.ll2010h 
and SDSS1035 ( South worth et al.ll2006 ') in quiescence. This lat- 
ter bright region was also detected i n Doppler maps of WZ Sge 
durin g its superoutburst in 2001 (Ba ba et ani2002t iHowell et alJ 
12003 ). We assume that the most probable interpretation of these 
tomograms is two spiral waves in the outer an nuli, formed by the 
2:1 re sonance in the extended accretion disk (iLin & Papaloizoul 
[T979I) . 

In January 2012, the intensity of the S-wave in the trailed 
spectra of the Ha line was somewhat fainter compared the total 
flux from the accretion disk. However, the distribution of matter 
in the accretion disk did not change overall. This observation is 
also confirmed by the map obtained at the same time (see 
FigurefTOt. 



7. Accretion disk structure of bounce-back systems 

Figure [TT| illustrates the current concept of CV evolution at the 
orbital period turn-around point and displays the position of the 
bounce-back candidates on the mass-transfer rate and mass-ratio 
to orbital period diagrams. Some CV classifications are are la- 
beled on the mass transfer rate to orbital period diagram. The 
classification in astronomy is often phenomenological: CVs with 
orbital periods close to the 80 min orbital period minimum that 
undergo infrequent (years to decades) superoutbursts are called 
WZ Sge-type stars. Objects with short periods that have not been 
observed in outburst or superoutburst but have spectral charac- 
teristics similar to WZ Sge (i.e, a mildly blue continuum with 
relatively weak Balmer emission lines from the accretion disk 
embedded in a broad absorptions formed in the atmosphere of 
the white dwarf) are listed as WZ Sge-type candidates. Bounce- 
back systems are spectroscopically similar to these, but not every 
WZ Sge-type object has necessarily passed through the turning 
point. The parameters of bounce-back systems (i.e. objects that 
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Fig. 12. Size of a typical bounce-back system and its compo- 
nents in solar units compared with the size of Jupiter and the 
sun. The 2:1 and 3:1 resonance radii are also shown. 
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Fig. 13. Primary's magnetosphere radius for different mass 
transfer rates (in solar mass per year) vs. the primary's magnetic 
field. 



would fall into the box in the right panel of Figure[TTJ are listed 
in Table E] 

To separate the bounce-back systems from the others, the 
mass of the secondary or the mass ratio is needed. Direct mea- 
surement of the mass of either stellar component in CVs is an 
extremely difficult task. But the mass ratio can be estimated 
from the period of the superhumps that may be seen in the light 
curves during superoutbursts. Four objects (V455 And, AL Com, 
SDSS0804 and EG Cnc) produced WZ Sge-type superoutbursts 
in the near past. Estimates of the system parameters often show 
the presence of a massive white dwarf (Mi > 0.9Mo) and an 
extremely low mass ratio q < 0.06, which assumes that the sec- 
ondary is a Jupiter-sized brown dwarf (see Figure [T2l i. A q this 
low implies that the Roche lobe of the primary is huge and the 
accretion disk could extend to the 2:1 resonance radius and be- 
yond. If that is the case, the disk will develop two spiral density 
waves in the outer annuli as shown by Kunze & Speith (2005]). 
The two waves will cause the appearance of the double hump- 
shape per orbital period in the light curves, especially in high- 
inclination systems. Thus, one way of separating bounce-back 
systems from other WZ Sge objects is to detect the spiral struc- 
ture in the accretion disk in quiescence. WZ Sge usually exhibits 
double-humps only during the early stages of the superoutburst, 
when the disk expands. In the left panel of Figure[TTl the mass 
ratio leading to the 2:1 resonance disk radius is indicated by a 
tiny dotted line, which clearly separates bounce-back systems 
from the rest. 

We aim to determine the conditions in bounce-back sys- 
tems that enable detection of these two spiral waves in the 
outer annuli. The accretion disk structure of bounce-back can- 
didates is analyzed here primarily based on ou r spectral an d IR 
photometric observations of SDSS1238 (.Aviles et al.ll2010h and 
SDSS0804 in quiescence. The conditions of the disk are inferred 
from the Doppler tomography maps and from the analysis of 
the SEDs of these systems. The SED modeling requires a low 
contribution from the accretion disk to the total flux. The main 
uncertainty of the SED model is the slope of the accretion disk 
continuum. The shape of a standard spectrum (Fx ~ /l"^^^) of the 
accretion disk is based on a blackbody approximation of a disk 
element flux intensity and the relation Teff(r) w r,(r/7?i)"^'''* for 
the radial tem perature structure of a steady-state accretion disk 
dWarnedll995l) . However, an optically thick accretion disk emit- 



ting as a blackbody of this size will totally dominate the radiation 
of the system. Because observations show that the white dwarf 
radiation dominates in the optical range and the secondary's ra- 
diation dominates the JHK bands in SDSS0804 and SP SS 1238, 
the "standard accretion disk model' ' (iFrank et alJl2002h does not 
apply to bounce-back systems. 

It is known that accretion disks in CVs develop optically 
thin oute r regions for m ass transfer rates below about 10"^ 
Mo yr ' (IWilliamslll9"80h . fTvlenda (1981 ) confirmed this and 
added that an increase i n the radius of the or-model disk 
(IShakura & Sunvaevlll973h always increases the role of the thin 
region. For a < 1 the thin part of the disk is rather cool (< 
6000K). Because a large part of the accretion disk is optical thin 
in continuum (i.e., in the non-LTE regime and for a < 0.1), th e 
disk temperature can drop b elow 5000K dPumont et al.lll99ib . 
ICannizzo & Wheeled (Il984l) studied the vertical structure of a 
steady-state, a-model thin-accretion disk for an accreting object 
of 1 Mq. They found that, for low accretion rates, the disk struc- 
ture is optically t hin. For 0.01 < a < 1, the solution of disk 
equations (l)-(3) (ICannizzo & Wheeleii Il984h can be double- 
valued with high- (~ 5000 K) and low- (~ 2000 K) temperature 
branches. For a > 0.1 a warm solution is possible in the inner 
region of the accretion disk, but disk annuli at larger radii will 
be in a cold state with T < 2000/r. Only the low-temperature 
solution exists for a a; 0.1. As o- decreases with temperature, the 
tendency to develop cold solutions in quiescence is enhanced. 
Until now, models of disk emission spectra for these cool disks 
have not been calculated, the nearest approach is a flat spectrum 
in the range 3000- lOOOOA for flie mo del of the cool (T= 5000 K, 
a - 0.03) accretion disk dldan et alJ l2008). The next important 
aspect regarding accretion disks of bounce-backers is the non- 
uniform structure inferred from the Doppler tomograms and op- 
tical light curves alike. As akeady mentioned, the non-uniform 
structure was interpreted as spiral arms in the accretion disks 
(Aviles et al. 2010, see Doppler maps discussion and references 
therein). The inner part of the accretion disk in a standard model 
is usually optical ly thick and forms a continuum of the disk 
spectrum (Canniz zo & Wheeleii 11984 and references therein). 
However, for WZ-Sge-type systems, various authors have sup- 
ported the idea that the inner part of the disk is void during 
quiescence (see lKuuUcers et all ( 1201 lh : lMennickent et alj d2001h 
and reference therein). This void is assumed to explain the long 
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Fig. 14. Top panel: model configuration used to calculate the light curves of bounce-back systems. The bottom panel shows the 
vertical thickness of the accretion disk (left) and the temperature distribution (right) of the model. 



recurrence time for superoutbursts and the absence of normal 
outbursts. The missing inn er part of the accretion disk ma y 
be caused by evaporation (Mever & Mever-Hofmeister 1994^ 
or/and by magnetic field of the primary WD ( Matthews et all 
l2006l) . According to the equation = (GMwd)~^M~7 j^? 
where jd = Bq^^d ~ ^^^^ Gauss cm^ (Kluzniak & Rapoaport' 
l2007l and references therein), the size of the magnetosphere will 
increase with the decrease of mass transfer rate (see Figure[T3Tl. 
A relatively faint (< 1 MGauss) magnetic field is enough to form 
a cavity in the inner part of the accretion disk. Magnetic fields of 
about 0.5 MG and stronger are detected in magnetic CVs from 
the white dwarf spin modulation observed in optical light or X- 
rays. V455 And is a proven intermediate polar. I n SDSS0804 a 
number of high-frequency periods were detected dPavlenko et al.l 
1201 1). but none was identified as a spin period so far. X-ray ob- 
servations are needed to probe the magnetic properties of the 
white dwarfs in bounce-back candidates to verify the interme- 
diate polar hypothesis. A typical system with parameters corre- 
sponding to a bounce-back system may look like the sketch in 
Figure[T4](top panel). The simulated disk is thin in the inner re- 
gions and thick along the spiral density waves, which are located 
in the outer annuli. 

8. Light curve simulation 

An eye-catching double-hump-shape per orbital period in qui- 
escent light curves is a common feature of objects proposed to 
be bounce-back candidates. We developed a geometric model 
of CV disks with two spiral density waves in the outer annuli 
of a thin accretion disk to model light curves of bounce-back 
systems. The top panel of Figure[T4] presents a snapshot of this 
model. The system is comprised of a primary white dwarf, a 



Table 4. Parameters used for V- band light curve modeling 



^ mod 


i 






^' 


e* 


<5* 


d 




n 


Ro 


Ro 








n 


Nl 


75-80 


0.2 


0.45 


0.7 


0.06 


35 


0.2 


N2 


80 


0.1-0.2 


0.45 


0.9 


0.08 


25 


0.25 


N3 


75 


0.2 


0.45 


0.7 


0.06 


35 


0.05 


N4+ 


80 


0.2 


0.45 


0.7 


0.06 


35 


0.05 



Ml = 0.9Mo, q = 0.05, T, = llOOOK, T2=2000K, 
Port = 51005 = 1.42/i, a = 0.64Ro, Rli = O.49R0, 
M = 1.6 X 10-" Malyr,p' = 0.01 

* — parameters of the disk and the spiral arm corresponding to 
the model of the two-armed spiral structure in the accretion disk 
of^lachisu^tah (2004) (see equations 3,4 and 5). 
+ — the model with Mi = 0.7Mo 

i is the inclination; M| is the mass of primary; q is the mass ratio 
Ti and T? are the temperatures of the primary and the secondary 
a is the system separation; R^i is the radius of the primary Roshe 
lobe r„ and tout are the inner and outer radii of the disk 



secondary brown dwarf star, a stream of accretion matter, and 
an accretion disk. The white d warf is a spher e defined by the 
mass-radii relation (2.83b) froml Warnej (11995V The secondary 
is assumed to fill its Roche lobe, and the R oche lo be shape is 
calculated directly using equation (2.2) from I Warned (Il995h for 
equipotential 'I>(Li). The temperature of the primary and the sec- 
ondary are defined as initial parameters. Although illumination 
of the secondary by the primary in these systems is low, this 
illumination is also included. A moderate temperature gradient 
T{r) ~ Tifin) X r^^l'^ between the inner and outer edges of the 
disk and T ~ T(r) x (1 + Ox z(r)) from the spiral density waves 



10 



Zharikov, S., et al.: The accretion disk in the post period-minimum CV SDSS0804 



AL Com 



V455 And 



-0.2 
-0.1 

0.1 
0.2 





' 1 ' 1 


' 1 ' . 








/ 








1 , 1 


1 , " 







0.5 



1 

SDSS0804 



1.5 



Model Nl 



-0.2 
-0.1 

0.1 
0.2 



6=0.2 





1 














' \i=8(f/ - 


1 









0.2 
0.1 

0.1 
0.2 




Phase 



0.5 



U U.J 


1 1 ^ 

i L.J 

SDSS1610 


L 


- ' 1 

a □ □ " E 

- . 1 


' 1 ' 1 
1 , 1 


1 

1 


0.5 


1 1.5 


2 


i=80 


Model W2 


6=0.25 




/ V-N / 1- =0.2^ 






/ '"\V/ X =0.15\ , 


r ft Y 


. r 


\ y^/t =0.1 " v 

1 , 1 


1 


0.5 


1 1.5 


2 


M = 0.7 Mo 


Model m 


6=0.05 


' 1 

A/ 

" . 1 


' 1 ' 1 

1 , 1 


1 

^'^=80°^ 

1 



1.5 



Fig. 15. Light curves of AL Com, V455 And, SDSS0804, and SDSS1610 folded with their orbital periods are shown in the upper 
four panels. Simulated light curves from our models are shown in the lower four panels (see Table|4ll. 



(see Figure[T4| was assumed. Here, r(r,„) is the temperature of 
the internal radius of the disk, which depends on the mass ra- 
tio (see equations 2.35 and 2.36 in Warner (1995)), and is a 
free parameter that defines the temperature gradient in the spi- 
ral density waves. The model takes into account the positions of 
the bright structures in the Doppler maps, the large size of ac- 
cretion disk, and the description of the spiral density waves in 
[Hachisu et al. (2004). Figure[T4l presents the geometry used in 
the model (top panel) and grayscale images present the height of 
the disk (bottom left) and the temperature distribution (bottom 
right) of the accretion disk. 

We calculated a variety of models using the typical aver- 
age parameters of bounce-back systems presented in Table [3] 
In Figure[T5]four representative models corresponding to the en- 
tries of TableH] are plotted along with the observed, averaged 
light curves of ALCom, SDSS0804, SDSS1610, and V455 And. 
The main difference between the models is the choice of pa- 
rameters, which influences the spiral density wave structures in 
the disk. In the second model (N2), the spiral density waves are 
hotter and thus the amplitude is higher than in the model Nl 
(Tabled. To show the dependence on system inclination from 
the orbital plane and the radius of inner disk in the models, we 
plot three curves for each model in the lower panels of FigurefTSl 
Obviously, the inclination varies in the Nl model and the inner 
disk radius varies in the N2 model. For models N3 and N4, a sig- 
nificantly smaller Q is used for the Miy/j = O.9M0, / - 75° and 



A^WD - 0.1 Mq, i - 80° cases, respectively. Clearly, the double- 
hump-shape light curve is easily reproduced by these models. 
From our models, we find that the shape and amplitude of the 
variability depends on the parameters chosen for the two spiral 
density waves, the disk temperature distribution, the size of the 
inner cavity, and the outer radius of the disk. We also find that the 
amplitude of the double-humped variability decreases with de- 
creasing temperature gradient in the spiral density waves and/or 
the decline of the system's inclination. We note that changing 
the primary's mass has no effect on the shape of the light curve. 
However, as noted previously , a low-mass primary will not al- 
low the disk to reach the 2:1 resonance radius, which results in 
the spiral structures in the outer annuli of the accretion disk. 



The number of free parameters of the light curve model 
is large. Using a minimization to fit each light curve pre- 
cisely would be a very complicated task and is not the aim of 
this study. The simple models used in this work validate previ- 
ously reported ideas on the disk structure in WZ Sge systems - 
a large, cool, mostly optically thin accretion disk containing two 
dense/hot spirals. We add that these models offer a reasonable 
explanation for the observed characteristics of bounce-back sys- 
tems. 
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9. Conclusion 

We here continued our study of SDSS J080434.20+5 10349.2 fo- 
cusing on the nature of double humps in the light curve and their 
relation to the structure of the accretion disk. 

The object exhibited a;0.07 mag variability with 42.48 min 
period, exactly half the orbital period of the system, during 2008- 
2009 when it almost retuned to quiescent level after the March 
2006 superoutburst. In September 2010, the system underwent 
yet another superoutburst and approached its quiescent level by 
the beginning of 2012. The light curve once again showed dou- 
ble humps, but with a significantly smaller (~0.01mag) ampli- 
tude. No long-term SDSS1238-like variability and mini-outburst 
phenomena were detected. 

The SED of SDSS0804 in quiescence was fitted by the three- 
component model, which requires a massive >O.7M0 white 
dwarf with a surface temperature of ~12000K, a late-type brown 
dwarf, and an accretion disk with a weak contribution to the total 
flux. 

H„ and H/? Doppler maps obtained in quiescence from time- 
resolved spectroscopy persistently show excess emission from 
two opposing sides of the accretion disk. This pattern is a com- 
mon characteristic that SDSS0804 shares with SDSS 1238 and 
other bounce-back candidates. 

We constructed a geometric model of a bounce-back sys- 
tem to reproduce the observed light curves. The double-humped 
light curves in quiescence can be generated by the model in 
which the accretion disk extends to the 2:1 resonance radius, 
is cool (~2500K) and contains two spiral density waves in the 
outer annuli of the disk. The spirals should be denser and hotter 
than the rest of the disk. The inner parts of the disk should be 
optically thin in the continuum or be totally void. A synthetic 
Doppler map constructed using this model of the accretion disk 
can also reproduce the observed maps, as was recently shown in 
lAviles et al.l (120101. Fig.4) . 
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